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Variations in the spatial configuration of the interstellar magnetic field (ISMF) near 
the Sun can be constrained by comparing the ISMF direction at the heliosphere found 
from the Interstellar Boundary Explorer spacecraft (IBEX) observations of a 'Ribbon' 
of energetic neutral atoms (ENAs), with the ISMF direction derived from optical po- 
larization data for stars within ~ 40 pc. Using interstellar polarization observations to- 
wards ~ 30 nearby stars within ~ 90° of the heliosphere nose, we find that the best fits 
to the polarization position angles are obtained for a magnetic pole directed towards 
ecliptic coordinates of A, j3 ~ 263°, 37° (or galactic coordinates of £, b ~ 38°, 23°), 
with uncertainties of ±35° based on the broad minimum of the best fits and the range 
of data quality. This magnetic pole is 33° from the magnetic pole that is defined by 
the center of the arc of the ENA Ribbon. The IBEX ENA ribbon is seen in sightlines 
that are perpendicular to the ISMF as it drapes over the heliosphere. The similarity of 
the polarization and Ribbon directions for the local ISMF suggest that the local field 
is coherent over scale sizes of tens of parsecs. The ISMF vector direction is nearly 
perpendicular to the flow of local ISM through the local standard of rest, supporting 
a possible local ISM origin related to an evolved expanding magnetized shell. The 
local ISMF direction is found to have a curious geometry with respect to the cosmic 
microwave background dipole moment. 

Subject headings: ISM: magnetic fields, clouds, HI — solar system: general — stars: 
winds, outflows 



Introduction 



The recent unexpected discovery by the Interstellar Boundary Explorer mission (IBEX. lMcComas et al 



2009b() of a 'Ribbon' of energetic neutral atoms (ENA) provides the first direct evidence, and in- 



ized cloud surrounding the Sun ( 


McComas et al. 


2009a; 


Fuselier et al. 


2009; 


Funsten et al. 


2009; 


Schwadron et al. 


2009; 


Mobius et al. 


2009b. The IBEX Ribbon is visible at 1 AU for directions 



where the sightline is perpendicular to the interstellar magnetic field draping over the heliosphere. 
The center of the Ribbon arc defines a value for the direction of the interstellar magnetic field 
(ISMF) in the interstellar material (ISM) that is interacting with the heliosphere. The ISMF di- 
rections traced by the IBEX Ribbon, and weak interstellar polarizations observations of towards 
nearby stars, together provide a unique constraint on the local ISMF direction, the outer boundary 
conditions of the heliosphere, and the global magnetic turbulence in the low density interstellar 
material near the Sun. In this paper we compare the available data on the ISMF direction from 
optical polarization of local starlight, with the ISMF direction signified by the IBEX Ribbon. 
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Interstellar material (ISM) sets the outer boundary conditions of the heliosphere, including the 
therm al and ram pressures of the gas, and magnetic pressures (|Davislll955l : [Parked 1 1 96 ll : iHolzer 
1989Q . The heliosphere nose is defined by the upwind direction of the flow vector of interstellar 
gas through the heliosphere, and is 15° above the galactic center and 5° above the ecliptic plane. 
The 0.05-6 keV (~ 50 — 1000 km s _1 ) ENAs detected by IBEX are formed in regions where 
charge-exchange between interstellar neutral hydrogen and energetic ions create inward flowing 
particles that reach the IBEX detectors. The ENAs formed inside the solar wind termination shock 
propagate away from the Sun, but the subsonic solar wind and pickup ions outside of the termi- 
nation shock, with momentum components dire cted back towards the Sun, can create the ENA s 
measured in the inner heliosphere by IBEX (e.g. Izmodenov et all2009 : Heerikhuisen et al. 2008 ). 
The subset of ENAs that propagate back towards IBEX at 1 AU are proving an excellent diagnostic 
of the heliosphere boundary conditions. The center of the Ribbon arc indicates an interstellar mag- 
netic field (ISMF) di rection towards ecliptic coordinates A,/3= 221°,39° (or in galactic coordinates 
towards ^,6=33°, 55°, Funsten et al.l 12009c Fuselier et al.ll2009|) . The Ribbon's arc does not follow 
the equator of the interstellar magnetic fie ld (ISMF) defined by th e arc center, but rather traces the 
distortion of the ISMF by the heliosphere (|Schwadron et al.ll2009l) . For magnetic field B and radial 
direction R, the length and width of the B ■ R region depend on magnetic pressure and the limit 
set on B ■ R, respectively. The Ribbon latitude, in comparison to the magnetic pole defined by the 
arc center, depends on the distance beyond the heliopause of the Ribbon formation. 

The ISM in the low column density sky is formed into filamentary anisotropic structures, 
whether it is neutral gas seen in the 21 -cm H° line, ionized gas observed in Ha emission, or dust 
traced by infrared emission or the optical polarizations of starlight. Magnetic fields are a candidate 
force for controlling this observed filamentary structure, since they generate anisotropic forces on 
the ISM through pressure that acts perpendicular to the field lines, and tension that acts parallel 
to the field lines (|Heilesl I20Q9Q . Prior to the heliospheric diagnostics of the very local ISMF, such 
as the IBEX Ribbon, the only data capable of constraining the ISMF within ~ 50 pc were high- 
sensitivity measurements of the polarization of starlight by magnetically aligned interstellar dust 
grains. The IBEX Ribbon is a feature that can be simulated with full heliosphere models, that 
simultaneously give the pressure of the ISM, and the orientation and strength of the ISMF. In the 
upwind direction (near the galactic cente r), the Sun appears to be at or close to the boundary of 



the s urrounding cloud (|Frisch et al.ll2009|) . The ENA production simulation of iHeerikhuisen et al 



(120101) . which is unconfirmed, predicts that IBEX is capable of detecting variations in the Ribbon 
configuration that would be caused by variations in the interstellar densit y or magnetic field , such 
as expected when the Sun exits the cloud now around the heliosphere (iFrisch et al.ll2010l) . The 
relation between the ISMF directions traced by the IBEX Ribbon, and the observations of the weak 
interstellar polarizations of light from nearby stars, provides a powerful tool for understanding both 
the distortion of the ISMF at the heliosphere, and deviations from a homogeneous ISMF field over 
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the nearest tens of parsecs. 



The magnetic configuration traced by the ENA ribbon, together with the asymmetry of the 
heliosphere, is a function of angle between the interstellar gas flow vector and the ISMF direction. 
The flow of interstellar He° through the heliosphere was measured by the Ulysses probe, and 
showed a warm cloud T ~ 6300 ± 340 K, with velocity 26.3 km s -1 and an upwind direction 
towards A,/3~ 255.5°, 5.1° (or £.6=3.6 o .15.1°. lwittell2004l after converting to J2000 coordinates). 
The Ribbon data therefore indicate an angle of ~ 46° between the He° flow velocity and field 
direction. The gaseous component of the heliosphere boundary conditions has been determined 
from models of the cloud ionization state based on interstellar absorption lines, which indicate 
that the surrounding inte rstellar gas is partially ionized with hydrogen ~ 25% ionized, and low 
density, n tot ~ 0.27 cm" 3 ( Slavin & Frischll2008 . e.g. model 26). The ratio of thermal and magnetic 
pressures, (3, is unknown for the surrounding cloud. If the local ISM is near the boundaries of an 
evolved superbubble, which is not yet clear, /3 ~ 1 might be expected. For the assumption of 
(3 = 1, the interstellar field strength is ~ 2.8 /iG (model 26). For comparison, Zeeman splitting of 
H° 21 -cm lines show that superbubble shells with filamentary structures have high ISMF st rengths 



and s mall p < 0.1, although column densities are over 100 times values in the local ISM (IHeiles 



1989b . Low (3 values are also found for the Local Bubble walls foreground of the Coalsack Nebula 
based on a Chandrasekhar -Fermi analysis of the polarization component in the plane of the sky 
(kndersson & PottejEooa) . 



The ISMF direction at the heliosphere is predicted by self-consistent heliosphere models 
(e.g. iRatkiewicz etaD boool : IPogorelov etaD l2004i 120091 : bpher et all 120091 : llzmodenovl l2009h 
that are constrained by the interstellar neutral and plasma density in the circumhelio spheric ISM 
(ISlavin & Frischl 120081) . the 10 AU difference between the Voyager 1 and Voyager 2 measure- 
ments of the term ination shock distance (of which 3-4 AU could arise from variable solar wind, 
Stone et al.ll2008 ), and the ~ 5° offset between the d irections of interstellar H° and He° flowing into 



the heliosphere (Witte 2004; Lallem ent et al 



a common coordinate system). 



2005 



Heerikhuisen et al 



, after precessing the He° and H° directions into 
(12010 ) models reproduce the IBEX Ribbon for 

223°, 40° (or l,b~ 35°, 54°, hereafter i,b are galactic co- 

30° ±5°, 33° ±4.3°, 



an ISMF direction directed towards X,(3 

ordinates). This direction is in good agreement with the ISMF direction £,b 

for B ~ 3.8 /xG, fro m MHD models of heliospheric asymmetries of IRatkiewicz & Grygorczuk 



(I2008|) . lOpher et alJ (I2009|) use the heliosphere asymmetries from Voyager data and the H°-He° 
offset to conclude that the ISMF field direction is inclined by 60° — 90° with respect to the galactic 
plane. Several heliosphere models suggest t he ver y local ISMF is directed towards the southern 
solar magnetic pole, such as iSwisdak et al.l (|2010|) who base their conclusion on the assumption 
that the 3 kHz emissions detected by the Voyager satellites during the early 1990's are formed by 
reconnection. The general agreeme nt between the ISMF determine d from the ENA ribbon, the 
Heerikhuisen et al. model, and the IRatkiewicz & Grygorczuk! (|2008l) model suggest directions for 
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the ISMF at the heliosphere near the center of the Ribbon arc. 



The new Ribbon diagnostic of the ISMF at the heliosphere provides an opportunity to evaluate 
the large-scale distortion or turbulence of the local ISMF, e.g. over ~ 40 pc, through comparisons 
with the polarization of starlight caused by magnetically aligned interstellar dust grains. Optical 
polarization data show that once the curvature of spi ral arms is ta ken into account, the global 
ISMF beyond ~ 0.3 kpc, is oriented towards £~ 83° (IHeileslll996l). Faraday rotat ion data show 
that the polarity of the global ISMF is directed towards £~ 83° (ITaylor et all 120091) . The Loop I 
magnetic superbubble, centered about 100 pc away, is a local large distortion of the global field 
with an angular diameter of ~ 160° (e.g. iHeilesI 1 1 976L 1 19981 : IWollebenl |2007[) . It is prominent 
in dust, gas, synchrotron emission, optical polarization, and Faraday rotation data. Loop I is the 
best candidate for the phe nomena that links the ISMF at the heliosphere with the ISMF causing 
nearby star polarizations (|Frischll2010i) . The overall geometry of Loop I, and the 'SI' subshell, 
indicate that S 1 has expanded to the solar location if it is approximately spherical, and the ISMF 
direction appears consistent with nearby polarization data for a field direction near 71°, 48°, 
with uncertainties of ~ ±30°. 

The ISM within ~ 20 pc of the Sun has very low column densities, typically < 10 18 5 cm~ 2 , 
and therefore low extinctions. Pioneering efforts to measure polarized starlight caused by mag- 
netically aligned interstellar dust grains found that nearby space is very empty of the ISM, except 
for a patch of nearby dust, within 40 pc, primarily in the fourth Galactic quadrant between 270° 
and 360°, and in the southern hemisphere (|Piirolall 19771 : lTinbergenl ll982). This dust 'patch' has a 
very local component, as the 2.5a detection of polarization towards the star 36 Oph, 5 pc beyond 
the heliosphere nose, shows. (It is remarkable that many of these weak polarizations line up along 
the IBEX Ribbon, see Fig. 1) The cloud giving rise to the polarization towards 36 O ph is not the 
circumhelio spheric ISM, since the velocities of the two clouds differs by ~ 2 km s _1 (IFrisch et al. 



20091) . so that the polarizations in this patch potentially trace a different ISMF direction t han the 
Ribb on arc. The next set of studies of the polarization of nearby starlight was conducted by LLeroy 
(|1993|) . from the northern hemisphere. Leroy made a catalog of observations of 1000 stars within 
50 pc, with different sensitivity levels, and reconfirmed the emptiness of nearby space. However, 
he was unable to either confirm or disprove the existence of the polarization patch observed by Tin- 
bergen. More recently, polarization observations with high me asurement accurac y have been made 
with the PlanetPol instrument, primarily in the northern sky (|Bailey et al.ll2010|) . The Tinbergen, 
Piirola, and PlanetPol data are used in this study. 

The goal of this study is to test the ISMF direction indicated by the center of the arc of the 
IBEX Ribbon with the ISMF direction traced by the polarization of light from nearby stars, using 
mainly polarization data from the literature with a range of sensitivity levels. These results are 
preliminary, in the sense that more and better high-sensitivity polarization data may affect the 
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conclusions. 



2. Local Magnetic Field Direction from Interstellar Polarizations 



Magnetically aligned dust grains in the ISM create a birefringent med ium with lower opacities 
parallel to the ISMF direction (e.g.. Davis & Bergdll968l : lLazarianll2003i) . The position angles of 
weakly polarized light, < 0.02%, from nearby stars can then be used to trace the local ISMF direc- 
tion. A distance limit of 40 pc is selected for polarization measurements of the local ISM, since the 
"boundary" o f the Local Bubble i s at about 50 pc for galactic longitudes 300°-360° and positive 
latitudes (e.g. lBarstow et al.lll997f) . We assemble polarization data for nearby stars, omitting stars 
with known intrinsic polarization or circumstellar disks, and systematically evaluate the polariza- 
tion position angles for a grid of i possible ISMF directions, ks,i, hi order to determine the ISMF 
direction (e.g. rotated coordinate system) that provides the best alignment between polarization 
position angles and the meridians of the ISMF. The polarization data used for these comparisons 
consist of three sets of archival starlight polarization data, and new unpublished observations for 
three stars. The star s a mple is res tricted to obje cts within ~ 40 pc of the Sun. The combined 
data sets of Tinbergen ( 1982 ^1 and Piirola ( 1977 ) provide data on ~ 140 stars, at la sensitivity 
levels of ~ 7 x 10~ 5 degree of polarization (or equivalently 0.007% or 70 parts per million), and 
in both northern and southern hemispheres. The strongest polarizations are seen towards the he- 
liosphere nose region, with st rengths of <y 0.02%. Recent PlanetPol measurements at sensitivities 



of a few parts per million by iBailey et al.1 (120101) add additional data in the northern galactic and 
ecliptic hemispheres, where polarizations are typically very weak (< 0.004%). Unpublished ob- 
servations of several stars are also available, including of (3 Oph (HD 161096), with polarization 
P = 0.00506 ± 0.00008%, and equatorial position angle fl RA = 166.2° ± 0.8° (S. Wiktorowicz; 
data acquired with the POLISH instrument, Wiktorowicz & Matthewsll2008 ). where 6> R a is the po- 
larization position angle in equatorial coordinates. We add measurements of polarizations towards 
A Sgr (HD 169916) and r Sgr (HD 177716) , that were acquir ed by A. Berdyugin and V. Piirola 
in the R-band with the DIPOL polarimeter (IPiirola et al.ll2005|) on the KVA-60 remotely operated 
telescope. For A Sgr, P = 0.036 ± 0.007% and 6 RA = 109° ± 7°. For r Sgr, P = 0.028 ± 0.007% 
and 0ra = 113° ±° 7. There are several stars with data from multiple sources (e.g. A Sgr, r Sgr, 
and (3 Oph). In these cases all of the data were included with equal weight in the analysis. 



The optical polarization position angles are plotted in galactic and ecliptic coordinates on 



'During 1973-1974 when the southern hemisphere Tinbergen data were acquired, the solar magnetic polarity was 
north pole positive (A>0, field lines emerging at the north pole). The solar polarity should have no effect on the optical 
polarization, as long as the grains are truly interstellar and outside of the heliopause. 
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Aitoff projections in Figs. [Tj [U The plotted size of the polarization vector is not related to the 
polarization strength, which spans over an order of magnitude for these stars (see Fig. [5]). The 
position angle of a vector is defined with respect to the north-south meridian passing through the 
position of the object, with position angle increasing in the direction of increasing longitude. 

To derive the orientation of the local ISMF we used a minimization procedure in which we 
assumed that the local ISMF has a dipole configuration, so that the variations in the orientations of 
the observed polarizations are due to the location of the field poles in the sky. If this assumption is 
correct, there should be a coordinate transformation that transforms all of the observed polarization 
vector dir ections into vectors that are parallel to a meridian of the transformed coordinate system. 



Following I Appenzellen (|1968[) . the position angle is transformed to coordinate system k B ,i using 
the relation: 

cos(k) * tan(b N ) - cos(^ - £ N ) * sin{b{) 

COt{6KA ~ 9i) = sinft - 4) (1) 

where 6> R a is the position angle in the equatorial coordinate system, and 6 h l h and b\ are the 
polarization position angle and star coordinates in the i th coordinate system corresponding to k B ^. 
The north pole of the equatorial coordinate system is located at Z N , 6 N in the i th coordinate system. 
Eq. 1 is used to transform between both equatorial and ecliptic or galactic coordinates, and between 
equatorial coordinates and the rotated ISMF frame k B ^. 

The ISMF direction that provides the best fit to the polarization position angles is selected 
by testing the polarization vector directions against the k Bi grid of possible ISMF directions. It 
is assumed that the "correct" ISMF direction will be parallel to the polarization vectors. The grid 
of possible ISMF directions are spaced by 1° in galactic longitude, £, and latitude, b, over the 
intervals 0° < l\ < 360° and < 85°. This comparison between polarization position angles 
and the ISMF direction was restricted to the subset of stars that have polarizations larger than 3a, 
where la is the measurement uncertainty, and that are located within ~ 90° of the heliosphere 
nose. Since the heliosphere nose is ~ 15° from the galactic center, the fitting region is essentially 
restricted to the first and fourth galactic quadrants where ISM and the magnetic sky are dominated 
by the magnetic shell of Loop I. Polarization position angles for each star j were first rotated into 
the k Bi coordinate system using eq. 1, to obtain the new position angle fly relative to the new 
ISMF direction corresponding to the grid point k B ,i- Polarization position angles are degenerate 
with respect to the north and south meridian directions, and the values of fly were corrected to be 
between 0° and 180° by adding or subtracting 180°. 

The best fitting ISMF then becomes the direction where the rotated ISMF (the i th grid point) 
yields a minimum of some function that describes a good match between the rotated polarization 
position angles and the rotated field direction. The function for determining the best fitting ISMF 
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that is adopted here is the direction where the ISMF pole corresponds to the i grid point that yields 
the minimum value for the function F h where = sin(9- h i) is the mean sin of j polarization 
position angles fly in the rotated coordinates. All data points are weighted equally, rather than by 
measurement accuracy, since otherwise data sampling the Ribbon region would have lower weights 
since the older data (from Tinbergen, 1982, see Fig. 1) have lower measurement accuracies. The 
position angle is the angle between the polarization vector and a meridian, so that for a perfect 
fit the rotated position angles will be 0° or 180°, and F{ = 0. Figs. 13141 shows the color-coded 
map of F\. The best fitting ISMF direction, k B best , that is determined where the minimum of 
Fi = 0.46, is directed towards the ecliptic coordinates A, (3 ~ 263°, 37°, or galactic coordinates 
£, b ~ 38°, 23° (Table 1). The broad minima shown in Figs. I3I4[ together with the range of data 
quality, suggest that a more accurate uncertainty for the best fit is ±35°. Figs. 13141 also displays 
(only) the polarization position angles that were used in the ISMF fitting process. 

Several tests of the fitting process were made. When the stellar data set was restricted to stars 
within 35 pc, the best fit ISMF direction changed by ~ 10° — 20° because four measurements near 
the ecliptic equator were removed from the sample, leaving a bimodal sample biased towards stars 
in the northern ecliptic hemisphere. The fitting function based on F-" = sin{6^) was also tried, 
and it gave a best-fitting ISMF direction towards £,b=53°, 25°, however this function overweights 
position angle values near the equator of the rotated ISMF, and so seems less suitable for this small 
set of data. The fitting procedure was repeated by varying the initial coordinate system used in the 
fit (e.g. equatorial, ecliptic, or galactic), and the results agreed to within ~ 1°. Another check was 
made by omitting the step of converting the rotated fly values to the range of 0° — 180°, yielding 
as expected the same result. 

The polarization position angles in the rotated frame have been constrained to be between 0° 
and 180°, so that the best fit value F min = 0.46 corresponds to mean position angles of 27°. When 
the standard deviation of the position angles is included, the best-fit mean position angle in the 
rotated frame is 27t.il degrees. In principle the dispersion in the position angles for the best-fitting 
ISMF direction could either be due to variations in the global configuration of the nearby ISMF 
over scale-lengths comparable to typical scales of energy injection, or to small scale turbulence at 
scale lengths typical of the plasma and magnetic properties of the partially ionized gas. However 
the intrinsic measurement accuracies of the data from the northern versus southern hemisphere 
data sets differ substantially, so that understanding small scale magnetic turbulence will require 
higher precision data in the southern hemisphere. 

For a uniform distribution of interstellar dust near the Sun, and constant grain alignment 
efficiency, the polarization strengths will increase as the angular distance between the star and 
the poles of the ISMF increase, i.e. the polarizations are strongest where the sightline is more 
perpendicular to the ISMF direction. These data do not show such an effect (Fig. [5]). Instead, 
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Table 1 . Best-fitting Magnetic Field Pole 



Coordinate 


Longitude 


Latitude 


System 






Polarization data 


- interstellar magnetic field: 




Galactic 2 


38° 


23° 


Ecliptic 


263° 


37° 


Center of Ribbon 


arc: 




Galactic 


33° 


55° 


Ecliptic 3 


221° 


39° 



Galactic coordinates are denoted by £, b , and ecliptic 
coordinates by A, (3. The estimated uncertainties on the best 
fit direction are ±35°, based on the broad minimum for the 
best-fit function, F;. 

2 This direction makes an angle of ~ 71° with respect to 
the vector motion of the flow of ambient local ISM past the 
Sun, in the LSR, which is from £,6=33 1°, -5° with velocity 
of -19.4 km s" 1 (iFrisch & Slavmhood) . 



3 This direction makes an angle of ~ 46° with respect to 
the heliocentric vector motion of the flow of interstellar He° 
into the heliosphere, w hich is from A,/3~ 255°, 5° with ve- 
locity of -26.3 km s" 1 Jwittell2004h . 
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stars with ecliptic latitudes j3> +10° consistently show much smaller polarizations than stars with 
ecliptic latitudes below +10°. This statement is also nearly true when galactic latitudes are used 
instead. All stars with polarizations less than 0.01% have ecliptic latitudes greater than /3=10°. 
All stars with polarizations larger than 0.01%, except for HD 150997, are located at more negative 
latitudes, (3< 10°. This effect follows from the distribution of ISM very close to the Sun, within ~ 
15 pc, which has higher column densities towards nega tive galactic latitud es than towards positive 



galactic latitudes in the galactic center hemisphere (e.g. lFrisch et al.1 12009). The ecliptic latitude of 
HD 150997 is +60°, and it is 26° from the ISMF pole at £=38°, 6=23°. A single isolated clump of 
dust towards HD 150997 is possible, or the polarization may be intrinsic to the stellar system. In 
Fig. [5] the stars used in the fit are color-coded according to the data source. Other significant data 
points, with polarizations larger than the 3a data uncertainties but not used, in the fit are plotted as 
open squares. 

Based on the above discussions we estimate uncertainties of ±35° on the best-fit ISMF direc- 
tion of A, j3 ~ 263°, 37° in ecliptic coordinates, or £, b ~ 38°, 23° in galactic coordinates. This 
direction is 33° from the ISMF direction at the heliosphere determined from the arc of the IBEX 
Ribbon. 



3. Discussion 



The ISMF direction of £, b ~ 38°, 23°, found from local polarization data, is directed towards 
the tangential region of Loop I at the North Polar Spur (NPS), that rises 50° vertically from the 
galactic p l ane ne ar £~ 30° — 40°. The NPS is a region of radio intense synchrotron emission, that 
Wollebenl (120071) attributes to the collision of the S 1 and S2 subshells of Loop I. Positive Faraday 
rotation measures for distant pul sars and extragala ctic sources towards the NPS indicate an ISMF 
that is directed towards the Sun jTaylor et al. 2009 ). However at the southern galactic latitudes for 
the same longitudes, where Wolleben's model suggests the SI subshell has expanded to the solar 
location, Faraday rotation measures are negative as is consistent with an ISMF d i rectio n pointing 
away from the Sun and towards the azimuthal field direction of £~ 83°. ISalvatil (|2010|) analyzed 
Faraday rotation measure and dispersion data towards four pulsars, 150-300 pc away in the low 
density third galactic quadrant, and found an ISMF directed towards £,b= 5°, 42°, with strength 
3.3 /iG, and with polarity directed into the northern galactic hemisphere. These data give the ISMF 
in the downwind direction, and it is within 33° of the best-fitting polarization ISMF direction and 
within 22° of the ISMF direction indicated by the center of the Ribbon arc. The good agreement 
between these three independent methods of obtaining the ISMF direction suggests the three kinds 
of measurements are tracing the same ISMF, and that it is relatively smooth. 



The formation mechanism for the Ribbon must be understood in order to relate the magnetic 
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field direction derived from the Ribbon configuration, with the magnetic field derived from po- 
larization data. The 26.3 km s 1 relative motion between the heliosphere and circumheliosphere 
ISM displaces the ISMF so that it drapes over the heliosphere, and the geometry of the outer helio 



sphere depends on the angle between the ISMF and interstellar flo w vectors (e.g. lRatkiewicz et al 



20001 : IPogorelov et al J 120091 : lOpher et al.ll2009l : llzmodenovl 120091) . The center of the Ribbon arc 
defines a magnetic field d irection that makes an angle of ~ 46° with the gas flow vector (Table 1). 
Heerikhuisen et al] (12010) reproduce the Ribbon figuration using three-dimensional MHD models 
of the heliosphere plasma, coupled to interstellar neutrals described by a kinetic distribution, and 
ions with a Lorentzian distribution. In this model, the ENAs originate upstream of the heliopause 
in the region where the ISMF angle varies with the distance beyond the heliopause. The model 
reproduces the location of the IBEX Ribbon quite well for the ISMF direction of \,/3= 224°, 41°, 
which is close to the center of the Ribbon arc. The details of the Heerikhuisen et al. model are 
not yet substantiated, because the pitch-angle distribution of the underlyi ng pickup ion ring-beam 
may scatter oyer tim escales shorter than the charge-exchange lifetime (IMcComas et al.l l2009al : 
Florinski et al.ll2010h . These same MHD models also show that the Ribbon moves towards the 
equator of the distant ISMF as the ENA origin pushes further upstream of the heliopause. Sev- 
eral alternative sce narios have been discus s ed for the Ribbon formation, including an orig in in the 
inner heliosheath (|McComas et al.ll2009al : ISchwadron et all 120091: McComas et all 1201 01) . Com- 
parisons between magnetic field directions derived from the Ribbon arc and polarization data will 
yield information on large-scale magnetic turbulence in the solar vicinity, once the formation of 
the Ribbon is fully understood. 

The ISM towards the star 36 Oph (5 pc away and 10° from the heliosphere nose) p rovides in- 
sights into the ISMF and gas forming the polarizations observed in the upwind direction. iTinbergen 
(119821) observed a polarization of ~ 0.02% towards 36 Oph, a strength that is unusually high com- 
pared to the mean polarizations found over long sightlines. A single interstellar cloud, the 'G' 
cloud, is present in front of both 36 Oph (5 pc) and the nearest star a Cen (1.3 pc), with a velocity 
that differs by 2 km s _1 from the ISM velocity inside of the heliosphere (e.g. LFrisch et al.ll2009i 
and references therein). This suggests that there is a single polarization screen within 1.3 pc of the 
Sun in the upwind direction, and allows the possibility that the ISMF direction derived from the 
polarization data is sampling a different magnetic field than the Ribbon arc. Heliosphere models 
predict that the configuration of the IBEX Ribbon is sensitive to 20° variations in the direction of 
the ISMF, so that a solar transition be tween the two field directions should be readily apparent in 
the configuration of the IBEX ribbon (IFrisch et al.ll2010|) . The Ribbon models are not yet proven, 
and the uncertainties on the ISMF direction obtained from the polarization data are large. Never- 
theless, it appears that the different ISMF directions obtained from the center of the Ribbon arc 
and starlight polarizations can provide information on the true large scale variations in the local 
ISMF, and ultimately on small-scale magnetic turbulence. 
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The local ISMF direction derived from polarization data can be used to test the possibility 
that the ISMF is perpendicular to the flow of ISM past the Sun. The interstellar cloud surrounding 
the heliosphere is part of a cluster of local interstellar cloudlets (CLIC) th at has a mean flow 
velocity directed away from the center of the Loop I magnetic superbubble (e.g. lFrisch et al. l2009h . 



For comparing the kinematical CLIC with the geometrically defined configuration of Loop I, the 
interstellar velocities are first converted into the local standard of rest (LSR). The standard LSR 
conversion generally assumed for converting radio velocities to the LSR is used here, e.g. a solar 
apex velocity of 19.5 kms -1 , 56°, 23°. The result is an LSR upwind direction for the CLIC 
that is towards £,&=331°,-5°, and an LSR upwind direction for the circumheliospheric cloud that 
is directed towards £,6=318°,0°. Both of these directions are toward the central regions of Loop 
I. A value of ~ 71° is found for the angle between the mean LSR flow velocity vector of the 
CLIC and the best-fitting ISMF direction from polarization data, £, b ~ 38°, 23°. The ISMF 
direction is therefore nearly perpendicular to the cloud motions, similar to what might be expected 
for the expansion of an evolved superbubble shell in pressure equilibrium with ambient ISM. For 
the ISM around the heliosphere, the angle between the ISMF direction from the Ribbon arc and 
the heliocentric He° flow vector is ~ 46°. The spatial regions sampled by starlight polarization 
and CLIC interstellar absorption features are similar, since stars within 30-40 pc of the Sun are 
utilized in both comparisons. Additional information on magnetic turbulence in the flow of ISM 
past the Sun might be obtained from the Chandrasekhar-Fermi method, by co mparing the velocity 



dispe rsion of clouds with the dispersion of polarization position angles (e.g. lAndersson & Potter 



2006). 



There may be some contribution to the polarizations from interstellar grains in the outer he- 
liosheath, which will tightly follow the ISMF that interacts with and is deflected around the helio- 
sphere. The magnetic field upstream of the heliopause filters out grains with large charge-to-mass 
ratios, Q/M, and gyroradii that a re smaller than th e characteristic lengths between the heliosphere 
bow shock and heliopause (e.g. Jsiavin et al. 2009|) . The radii of the magnetically exclude d grains 



are < 0.01 — 0.1 fim, and are comparable to the sizes of interstellar polarizing grains (e.g. lMathis 



1986Q . The strongest polarizations occur for stars located on the northern edge of the IBEX Ribbon, 
which is 15° — 30° above the equator of the ISMF direction towards the arc center, corresponding 
to the blue points in Fig. [5J For truly interstellar polarization and a uniform dust distribution, the 
strongest polarizations are expected at the equator of the ISMF, so perhaps additional contributions 
to starlight polarizations from nano-sized grains in the heliosheath regions are possible. 

The lo cal ISMF dir ection found here has a curious coincidence with the CMB dipole (as found 



previously. iFrischl 120071) . The great circle midway between hot and cold poles of the cosmic mi- 
crowave background dipole moment bifurcates the heliosphere nose region and is aligned with the 
direction of the local ISMF direction found here, to within the uncertainties. The ISMF direction 

of i, b ~ 38°, 23° is at an angle of 90° ± 8° from the dipoles at £,b= 264°, +48° and 83°, -48°. 
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4. Conclusions 

In this exploratory study we find the best-fit to the polarization position angles towards ~ 30 
stars within 40 pc of the Sun and 90° of the heliosphere nose (16° above the galactic center), using 
mainly data from the literature, is £, b ~ 38°, 23° (or A, f3 ~ 263°, 37°, Table 1), with uncertainties 
of ±35° based on the flat minimum of the best fit direction. The ISMF direction indicated by the 
center of the IBEX Ribbon arc is £, b ~ 33°, 55°. The difference between the two, ~ 33°, is 
marginally significant given the uncertainties. The sensitivity of the IBEX Ribbon to variations 
of < 20° in the ISMF direction, and the fact that the Sun is in a cloud that is different from 
the nearest upwind ISM, support the possibility that the different ISMF directions obtained from 
the Ribbon and polarization data are tracking either large-scale distortion of the magnetic field 
direction, or possibly small-scale magnetic turbulence. This comparison is possible because the 
Ribbon is observed where the ISMF draping over the heliosphere is perpendicular to the sightline. 
Better constraints on the distortion of the ISMF in the solar vicinity will be possible once the 
formation of the ENA Ribbon is understood, and more high sensitivity polarization data towards 
the Ribbon are available. 

There are several implications of the ISMF direction found here. The ISMF vector direction 
is perpendicular to the bulk LSR velocity of the cluster of local interstellar clouds flowing past 
the Sun, which is consistent with an origin in an evolved expanding magnetized superbubble shell. 
The similarity of the ISMF directions found from the ENA Ribbon and optical polarization data 
suggests that the nearby ISMF is coherent over scale sizes of decades in parsecs, and that variations 
in the ISMF direction due to large scale distortion or magnetic turbulence are on the order of 35° 
or smaller. A curious coincidence is that the direction of the local ISMF is within 8° of the great 
circle that divides the hot and cold poles of the cosmic microwave background dipole moment, 
which also passes through the heliosphere nose. 
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Ecliptic Coordinates 
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Fig. 1. — The polarization vectors of stars within ~ 40 pc are show n in the ec l iptic c oordinate 
system, and c olor-coded for the data source. The ' 1974' data are from lTinbergenl (|1982L collected 
in 1974) and IPiirolal (|1977l) . The plot is centered on the heliosphere nose located at ecliptic co- 
ordinates (purple triangle) of A=255.4°, ,5=5.1°, and longitude increases towards the left in each 
figure. Symbol sizes do not indicate the strength of the polarization. The Compton-Getting cor- 
rected ENA fluxes at 1 keV are plotted for directions where the ENA count rates are larger than 
113 counts cm -2 s _1 sr - 1 keV" 1 , which is 1.5 ti mes the mean ENA flux at 1 keV as measured by 



the IBEX-HI instrument (McComas et al 



the center of the Ribbon arc (Funsten et al. 



2009a). The dotted lines show the ISMF determined by 
20091) . 
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Galactic Coordinates 
ISMF: L,B = 33°,55° 
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Fig. 2. — Same as Fig. [T] except the figure is in galactic coordinates, and centered on the galactic 
center. 
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Fig. 3. — The value of the function F\ = sin{6 h \) evaluated over a regular grid of % possible 
interstellar field directions (see $2]). The function is color-coded and plotted in the ecliptic coordi- 
nate system, centered on the ecliptic nose, at A ~ 255°. The gray dotted grid shows the best-fit in 
the ecliptic coordinate system to the ISMF, 5 i= best, which is directed towards A, (3 = 263°, 37°. 
Polarization vectors for the stars used in the fit are shown as either black or white bars, for visual 
clarity. 
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Fig. 4. — Same as Fig. [3l except quantities are plotted in galactic coordinates and centered on 
the galactic center. The best fit to the ISMF in the galactic coordinate system, B i=hcst , is directed 
towards £, b = 38°, 23°. 
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Fig. 5. — The polarizations of stars in the combined sample (in units of percentage of degree of 
polarization, ordinate) plotted against the angular distance between the star and the north pole of 
the best-fitting interstellar magnetic field, which is £~ 38°, 6~ 23°. The stars whose position 
angles are used to calculate the magnetic field direction are plotted as dots or triangles, and color 
coded to indicate the origin of the data. The stars in the full sample that have polarizations larger 
than the 3a uncertainties, but that were not used in the position angle fit (see text), are plotted 
as open squares. All stars with polarizations less than 0.01% have ecliptic latitudes greater than 
(3=10°, while all stars with polarizations larger than 0.01%, except for HD 150997, are located at 
more negative latitudes, (3< 10°. The ecliptic latitude of HD 150997 is +60°, and it is 26° from the 
ISMF pole at £=38°, 6=23°. 



